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Abstract 

Using FOPI data on elliptic flow in Au+Au collisions between 0.4 and 1.5A GeV we extract 
constraints for the equation of state (EOS) of compressed symmetric nuclear matter using 
the transport code IQMD by introducing an observable describing the evolution of the 
size of the elliptic flow as a function of rapidity. This observable is sensitive to the nuclear 
EOS and a robust tool to constrain the compressibility of nuclear matter up to 2 pq. 
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1. Introduction 

The equation of state (EOS) of nuclear matter is of fundamental interest and has been 
the object of intense theoretical efforts since several decades. The interest is boosted by the 
fact that it is an important ingredient in modeling fascinating astrophysical phenomena 
such as compact stars jl[ and core collapse supernovae j^. The calculation of the nuclear 
EOS from first principles, such as very recently attempted in Q, is a very complex task. 

The theory of nuclear forces, such as reviewed in j^, requires a confrontation with 
empirical facts, when it is applied to the ’ab initio’ problem of EOS calculations. In the 
last decade astrophysicists have pointed out (for a recent review, see [l|) that the EOS 
relevant for compact stars can in principle be deduced from the increasingly accurate 
and voluminous systematics of ’neutron’ star masses and radii by use of the relativis¬ 
tic Tolman-Volkov-Oppenheimer (TOV) ecmations. However, at the present time precise 
model-independent radii are still missing |^, while on the mass front line considerable 
progress has been made 0,0- Nevertheless the need of additional information remains: 
the composition of matter in the core of neutron stars is not known so far. 

A complementary method to approach the nuclear equation of state, practiced since the 
mid-eighties, is the use of heavy ion collisions over a wide range of incident energies, system 
sizes and compositions. The general difference between astrophysical and earth bound 
experiments is that the latter ones deal with rather isospin symmetric {N k, Z, N neutron, 
Z proton number) system compositions than are expected in the center of compact stars, 
and also, presumably, the temperatures in such reactions are much higher. There are, 
besides astrophysics, very good reasons to study symmetric (or almost symmetric) nuclear 
matter. It is clear that we must understand symmetric nuclear matter if we want to 
explain the properties of nuclei. Convincing theories are expected to reproduce properties 
of symmetric and ’pure’ (or ’almost pure’) neutron matter on the same footing. Also the 
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possibility of a phase transition as the density/temperatnre is raised can only be stndied 
qnantitatively in heavy ion reactions. 

However, deriving constraints on the nnclear eqnation of state (EOS) from heavy-ion 
collision data is a highly non-trivial task: we are in a dynamic sitnation with the character¬ 
istics describing the state of the system varying rapidly in time, while the observables are 
measnred long after the varions incoming and created parts of the system have separated 
forming a three-dimensional momentnm pattern hlling large (bnt not all) parts of phase 
space. Clearly, the claim to a signihcant ’constraint’ on the EOS cannot easily be made 
’mo del-independent ’. 

So far attempts to derive reasonably narrow constraints are limited to the energy regime 
below lOH GeV. The reason is that exclusively for this energy regime some kind of a clock 
is available. As described in a short, but pedagogical description [sj of the reaction scenario 
in this regime, a clock is provided by the fact that the ’spectator’ nucleons, not initially 
hit in collisions with non-zero impact parameters, are in contact with the more central 
’participant’ matter that forms the compressed ’hreball’ for a dehned time, the passing 
time, or “blockage time” j^. This time turns out to be comparable to the expansion time. 
The latter is governed by a fundamental property, the sound velocity which connects to 
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the EOS. 

The price one has to pay for a ’clock’ consists in more complex shapes 0, B 
due to the presence of spectators. In contrast, at energies beyond 10A GeV one can 
essentially start from ’simple’ almond or elliptic shaped participant clusters of nucleons 
whose evolution generating dominantly in plane ’elliptic’ flow was considered by Ollitrault 
0 in the early nineties. However the spectator clock is no longer operative. 

But even in this (nowadays ’low to intermediate’) energy regime the present knowledge 
on the properties of compressed nuclear matter is not satisfactory. 

Hopes that pion yields, measured in the mid-eighties, would be a textbook signal of 
a stiff EOS 12j have not materialized in the sequel. For possibly related reasons, the 


interpretation of pion isospin ratio data in helping to constrain the EOS of asymmetric 
systems is not settled [l^, . 

The use of positively charged kaons, , in the subthreshold region proposed by 
has been more successful. It is based on the principle that are produced in larger 
amounts if the nuclear EOS is soft due to the higher achieved density causing more fre¬ 
quent collisions. This effect is enhanced in the deeply subthreshold energy regime because 
the production of then requires multi-step reactions, including A production, in the 
medium. In addition, the relatively low reabsorption cross sections are supposed to favor 
a good memory of the high density phase. 

Using data from the KAOS Gollaboration on heavy ion systems of varying size (namely 


G-fG and Au-|-Au) and hence achieved density [ig, two theoretical groups [17|, [iSlj con¬ 
cluded that the so-called ’soft’ EOS (see later) was best suited to explain the data. Despite 
this encouraging success some questions remain open. The kaon is very rare at the lowest 
measured incident energy, 0.8A GeV, where the sensitivity to the nuclear EOS is largest 
and one has to ensure that all ’bulk’ observables (multiplicities, clusterization, stopping, 
flow) are under control. On the experimental side some desirable data, needed to complete 
the picture, are still missing: when switching from carbon to gold one is varying, besides 
the size, also the proton-neutron asymmetry, but the isospin partner, has not been 
measured at 0.8A GeV. Also, detailed experimental data on the nucleonic observables in 
G-l-G, needed to corroborate the simulated densities and momentum dependencies, are not 
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available. 

Clearly other observables such as particle flow are welcome in order to confirm the 
conclusions from the ’kaon method’. 

In a study [l^ of Au+Au collisions covering the beam energy range of 0.25A to 1.15A 
GeV the EoS Collaboration made an attempt to establish a ’complete characterization 
of squeeze-out’ (preferential emission out of the reaction plane) including a simultane¬ 
ous consideration of ’three categories of collective motion’ (radial and sideward flow and 
squeeze-out). From the work which was focused on the 0.6A GeV data, the authors con¬ 
cluded at the time (1996) that ’the QMD (a transport code j^) comparisons do not show 
a level of agreement with the data that justifies inferences about nuclear incompressibility’. 

Almost 10 years later (2005) the FOPI Collaboration published an excitation 
function of ’elliptic flow’ (see later) extending from 0.09A to 1.5A GeV beam energy for 
Z = 1 particles, also in the Au-|-Au system. The data were compared to the predictions 
from four different transport codes. The conclusion was ’no strong constraint on the EOS 
can be derived at this stage’. 

A transport theoretical study of sideflow and elliptic flow data (restricted to emitted 
protons) both from the BEVALAC and the AGS accelerators was undertaken in j^. The 
authors characterized their ’trial’ EOS by the incompressibility Kq of of symmetric nuclear 
matter around saturation density (p = po ~ 0.16 fm“^). An EOS labeled by Kq = 167 
MeV was adjusted to the experimental data on directed flow at incident energies below 
lA Gev, with a better description of the higher energy data using Kq = 200 MeV. The 
elliptic flow data seemed to require Kq = 300 MeV. 

In the present work we intend to improve the situation in the lA GeV regime. Profltting 


from extensive flow data published recently by the FOPI Gollaboration j22| we take here 


specifically a close look at the elliptic flow data, taking advantage of the afore-mentioned 
clock. In doing so 

1. we use not only protons, but also two- and three-nucleon clusters that are emitted 
in Au-|-Au collisions in the incident beam energy range from 0.4A to 1.5A GeV and that 
have larger flow signals than single nucleons 

2. we use not only mid-rapidity data, but use the information from 80% of the target- 
projectile rapidity gap; as we will show this information can be lumped into a single 
observable. 

The data are confronted to a well established transport code. Isospin Quantum Molec¬ 
ular Dynamics, IQMD j^, 24| using various phenomenological EOS to assess the EOS that 
is most compatible with the FOPI data. Four Skyrme type parameterizations are available 
in the code: H (’stiff’). S'(’soft’), HM (’stiff momentum dependent’), SM (’soft momen¬ 
tum dependent’). As in we characterize the EOS (for details see j^) by saturation 
density incompressibilities: Kq = 200 (’soft’) and Kq = 380 (’stiff’). 

We also estimate an error band from uncertainties in the data and compare to some 
recent representative microscopic EOS calculations. 

In the following section we present the analysis of the elliptic flow data obtained by 
FOPI and the resulting FOPI-IQMD constraint on the EOS. Dynamical details from 
the IQMD simulations are presented in section |3]and the constraints for the EOS obtained 
from the FOPI to IQMD comparisons are confronted with representative theoretical cal¬ 
culations in section m Finally, in section [5] we briefly discuss the obtained results and the 
needs for future efforts. 
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2. Analysis and results 


A comprehensive documentation on flow data in the SIS energy range measured by the 
FOPI Collaboration is available in j^. We shall therefore be brief here. 

Owing to collective flow phenomena, discovered experimentally in 1984 25, 0 , it is 
possible to reconstruct the reaction plane event-by-event and hence to study azimuthal 
correlations relative to that plane. The well established parameterization 27|, l28| of the 
observed azimuthal, (j), distributions is used: 


dN 
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Here we introduced the longitudinal (beam axis) rapidity y in the (c.o.m.) reference 
system and the transverse (spatial) component t of the four-velocity u, given by Ut = Ptl- 
The 3-vector /3 is the velocity in units of the light velocity and 7 = 1/ y/l — (3“^. Throughout 
we use scaled units yo = y/yp and Uto = Ut/up, with Up = Pp'jp, the index p referring to the 
incident projectile in the c.o.m.. In these units the initial target-projectile rapidity gap 
always extends from yo = to ?/o = 1- In 2^ the transverse momentum method |29t] 
was used including all particles identified outside the mid-rapidity interval |yo| < 0.3 and 
excluding identified pions. 

In Fig. [T]we show a sample of proton ’elliptic’, —V 2 , flow data from the FOPI Collab¬ 
oration (black dots with error bars) together with simulations using IQMD with a stiff 
version of the EOS (HM, red) and a soft version (SM, blue), see Fig. [3] for details to which 
we shall come back later. The reaction is Au-|-Au at an incident energy of 1 . 2 A GeV. In 
22 l | many more such data spanning the incident beam energy range 0.15 — 1.5A GeV and 


varying centrality and system composition and size are shown. 

By plotting —V 2 (rather than +V 2 ) one sees in this energy regime an enhancement of 
strength around mid-rapidity, i.e. predominantly out-of-plane emission, which was given 
the suggestive name ’squeeze-out’ in pioneering, and predictive, theoretical work 13 on 
the subject. Nowadays this phenomenon is named elliptic flow avoiding to conclude on the 
importance of shadowing and density gradients. At these energies the elliptic flow changes 
sign at high |yo| becoming predominantly in-plane. 

As in earlier FOPI work, collision centrality selection was obtained by binning distribu¬ 
tions of the ratio of total transverse to longitudinal kinetic energies in the center-of-mass 
system, ERAT. We estimate the impact parameter b from the measured differential cross 
sections for the ERAT using a geometrical sharp-cut approximation. We wish to stress 
here that we have used the same procedure (ERAT) for the simulations: the impact param¬ 
eter, not accessible directly by experiment, does not enter the analysis of the theoretical 
events. A similar procedure was used to correct for finite resolution of the reaction plane 
orient at ion ( see 22 j). 

We characterize the centrality by the interval of the scaled ’impact parameter’ 60 defined 
by 60 = b/bjnax-, taking b^ax = 1.15(Ap^ -|- A^^) fm. The centrality chosen for the figure 
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Figure 1: Proton elliptic flow data, —V 2 {yo), and IQMD-SM/HM simulations. See text for further expla¬ 
nations. 


is 0.25 < bo < 0.45. The choice of an optimum centrality for this study was guided by the 
intention to compromise between the need for sufficiently sized ’shadow bars’, composed 
by the spectator nucleons, on one side and the requirement of sufficient compression in 
the participant volume, a condition demanding a maximum size (maximum centrality) of 
the participant matter. 

Taking a closer look at —^ 2 ( 2 / 0 ) we see that the predicted shape is sensitive to the 
EoS in the full rapidity range. To take advantage of this feature we introduce a quantity 
dubbed V 2 n dehned by V 2 n = |'y 2 o| + |'y 22 | where the parameters are hxed by a £t to the 
flow data using ^ 2 ( 2 / 0 ) = ^20 + "^22 • 2/o scaled rapidity range I 2 / 0 I < 0.8. Such fits 

to the simulation data (here the error bars are of statistical origin) are shown in the left 
panel of Fig. [1] and then confronted with the experimental data (the error bars are now 
dominated by systematic errors) in the right panel. While V 20 alone characterizes the flow 
at mid-rapidity only, and the algebraic addition V 20 + V 22 determines the flow at 2/0 = 1 the 
quantity V 2 n combines features at both rapidities. The centrality is for 0.25 < Bq < 0.45 
and low momenta are cut off both in the data and the simulations. The cutting out of low 
transverse momenta, or equivalently, of scaled transverse 4-velocities {uto < 0.4), originally 
forced by apparatus limitations, actually turns out to raise the sensitivity of V 2 to the EOS 
as flow is generally converging to zero at low momenta (see e.g. 0 )- 

The V 2 n obtained for Au-|-Au between 0.4A and 1.5A GeV are shown in Fig. [2] for 
protons (lower left) and deuterons (lower right), tritons (upper left) and ^He. As the 
beam energy dependencies are rather weak, we indicate the average behavior by straight 
lines. The comparison of the data for V 2 n with the calculations shows a rather convincing 
preference for SM. The sensitivity is large: there is a factor 1.63 ± 0.06 between HM 
and SM, a difference exceeding significantly the indicated experimental error bars. If 
we compare this factor with the fluctuations of the experimental data points around the 
average values (dashed lines) we can estimate the uncertainty of the deduced EOS. 

To obtain quantitative conclusions on the EOS characterizer Kq for the energy range 
investigated in the present work we have constructed Table I which contains the relevant 
numbers (read off from the figure) with error bars. 

The Table lists the average (straight lines in the figure) V 2 n obtained, in order, for 
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Figure 2: Elliptic flow V 2 n for protons, deuterons, tritons, ^He as function of incident beam energy. 


SM, HM and FOPI, the ratio between HM and SM, the difference between HM and 
SM, the difference between FOPI and SM, and in the last colnmn the modification of K 
(for SM) needed to acconnt for the difference FOPI-SM, assnming for simplicity a linear 
interpolation between the Ko for HM and SM. From the information nsing fonr different 
particles we thns obtain a weighted average Kq = 190 ± 30, while the theoretical ratio 
between HM and SM is 1.63 ± 0.06. The nsed procednre is somewhat naive, bnt this is 
deemed nncritical considering the adopted error on Kq. 

We have plotted in Fig. Elthe resnlting EOS with its nncertainty band. For comparison 
we have also plotted the ’trial’ EOS, HM and SM, nsed in the simnlations. For the valid 
density range, indicated by vertical bars, see Section [3l Note the phenomenological EOS 
HM and SM inclnde the satnration point at p/po = 1, E/A = —16 MeV by constrnction. 
This hxes the absolnte position of the curves: the heavy ion data are only sensitve to the 
shape, i.e. the pressure which is essentially the derivative. Therefore the uncertainty of 
this fundamental point, about 0.5 MeV for the binding energy and 0.1 fm for pq, is not 
included in the uncertainty band. This will be important later when comparing to some 
representative theoretical EOS, section 01 However, we can already at this stage conclude, 
in complete agreement with ref. that a stiff EOS, characterized by Kq = 380 MeV is 
not in agreement with flow data in the incident energy range (0.4 — 1.5)A GeV. For OAA 


GeV this had also been suggested in [31 
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Table 1: Derivation of the effective K using the elliptic flow observable V 2 n 



V2n SM 

V2n HM 

V2n FOPI 

HM/SM 

HM-SM 

FOPI-SM 

AK MeV 


0.1658 ±0.0030 

0.2609 ± 0.0027 

0.1774 ±0.0028 

1.57 ±0.07 

0.0951 

±0.0116 

±22 ±8 

2 H 

0.3676 ± 0.0080 

0.6274 ±0.0087 

0.3237 ±0.0029 

1.71 ±0.08 

0.2598 

-0.0439 

-30 ±8 

3H 

0.5740 ± 0.0214 

0.9161 ±0.0252 

0.5223 ± 0.0048 

1.60 ±0.08 

0.3421 

-0.0517 

-27± 17 

3He 

0.5540 ±0.0217 

0.9048 ± 0.0265 

0.4537 ±0.0050 

1.63 ±0.16 

0.3501 

-0.1010 

-52 ± 18 


3. Simulations: the scenario 


In order to determine the EOS empirically, it is reasonable to start from low den¬ 
sities/pressures/temperatures and gradually raise these characteristics using sufficiently 
close steps. This is necessary to become aware of discontinuities (even modest ones) pos¬ 
sibly indicating important structural changes. 

The hrst step is to constrain the relevant density range. Isoscalar giant monopole 
resonance (GMR) data (’breathing modes’) have been used 3^ with the claim to determine 
the incompressibility around saturation density (p = Po = 0.16 fm ~^) where, by definition, 
the pressure in infinite symmetric nuclear matter is zero. An incompressibility K{pq) = Kq 


of 210 ± 30 MeV was deduced in [3^ from the monopole energies. 


In this context, it is interesting that recently some new insights have come up con¬ 
cerning the Kq value. In the authors have pointed out that the GMR actually tests 
compressibility at about 2/3po which is the average density of large nuclei whose hnite 
depth surface area is not negligible. The ’microscopic’ theories using energy density func¬ 
tionals (EDF), adjusted to reproduce the GMR data, then actually have to extrapolate 
to po leading according to |33| to Kq = 230 ± 40 MeV, the error reflecting fluctuations 
between different EDF. 

One should not expect that the value for the incompressibility of nuclear matter is 
directly comparable to the one used to characterize the tentative EOS adjusted to repro¬ 
duce heavy ion flows with transport codes, since the latter, as we will show, refer to data 
associated to densities significantly larger than po, and there is no basic principle which 
requires that the EOS is of simple parabolic shape over large ranges of density. 

Whereas the determination of the nuclear matter EOS with GMR data is an approach 
to saturation density ’from below’, our extensive IQMD simulations of the observations 
discussed in the present work suggest that in the FOPI experiments, we are advancing 
towards saturation density ’from above’. 

In contrast to model calculations done in an effort to reproduce the data, the simu¬ 
lations discussed in this section are done at a single impact parameter, since our goal is 
to elucidate the reaction scenario and to characterize which ’typical’ densities are probed 
with the proton elliptic flow that has been extracted from the FOPI experiments. For this 
purpose one has to determine at which times - in the course of the collision - and medium 
conditions are linked with and influence the most the development of the elliptic flow. 
Here, we restrict ourselves to the study of the ^®^Au-|-^®'^Au system at an ’intermediate’ 
impact parameter b = 3/m. It is close to the upper limit of the centrality range chosen 
for the experimental data in Figs. [T]and|2j 

In Fig. m various quantities are plotted as a function of time scaled to the passing time 
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HM/SM/FOPI 



Figure 3: Various symmetric nuclear matter EOS. Dashed (red) HM, full (blue) curve, SM. The dotted 
curves and the thin vertical bars delimit the FOPI EOS (crosses) constraints. The thick vertical bar is 
the density relevant for the GMR according to [33|. 

t/tpass, the time the nuclei pass each other is approximated by tpass = ‘^R/up, where Up was 
introduced in section [2J The time zero is dehned as the the moment projectile and target 
nucleons have the first hard collision. Typical values of tpass for various incident energies 
are listed in Table [2l The top panel of Fig. 0] shows on the left hand side the normalized 
density, p/po (where po is the nuclear matter ground state density) seen by the protons in 
the central region of the collision. The central region is dehned as a sphere with a radius 
of 1 fm around the center. This data is shown for reference and indicates the maximum 
density reached in those reactions. For a soft EOS higher densities are reached than for 
a hard EOS because the acting potential is less repulsive. The upper right panel depicts 
the mean density seen by all protons with transverse 4-velocities Uto > 0.4 and ||/o| < 0.8 
in the hnal stage of the reaction. This cut was adopted in the analysis of the FOPI data 
and we note that this selection removes most of the spectator particles. One observes 
that the maximum of density probed by these protons in the course of the collision is 
reached at around the full overlap of the incoming system {t/tpass = 0.5), and that the 
’soft’ EOS takes a slightly longer time to reach the maximum density, and reaches higher 







Figure 4: IQMD predictions of the time evolution of various quantities in Au+Au collisions at 1.5A GeV 
incident energy and an impact parameter 6 = 3 fm for various parameterizations of the EOS: soft EOS 
(SM) and hard EOS (HM), distinguishable through different line styles. The various quantities from left 
top to right bottom panels are: the reduced density probed by protons in the center of the colliding system, 
idem for all protons finally in the phase space region defined by |yo| < 0.8 and Uto > 0.4, the integrated 
(since time zero) sum of collisions experienced by these protons and their integrated momentum change 
due to the mean field (force). 


densities, as regard to the ’hard’ EOS. The mean time-integrated number of collisions 
for selected protons is shown in the lower left panel, and in the right panel their mean 
integrated momentum transfer caused by the potential. Obviously nearly all collisions 
have taken place until the passing time and most of the momentum transfer due to the 
potential has been achieved, i.e. the influence of the potential becomes very small because 
the system is getting diluted. We only see a slight dependence on the stiffness of the EOS 
in the number of collisions, caused by the different mean free paths of nucleons due to the 
different densities reached. However, the momentum transfer caused by the potential in 
the ’hard’ case is much larger than in the ’soft’ one, which reflects of the higher repulsion of 
the hard equation of state. Similar conclusions can be drawn for all incident energies that 
we have simulated (0.15 to 1.5 A.GeV). A summary of the maximum densities reached 
during the time evolution of the reaction for a range of beam energies is shown in Fig. [HI 
Fig. [5] synthesizes the dependence of the variable V 2 n as a function of reduced time 
t/tpass- This observable saturates only after two to three times the passing time with 
identical dependence for soft and hard EOS. To underline the similarity we applied a 
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Figure 5: IQMD results for V 2 n (see Section as a function of scaled time t/tpass for different EOS. 
The values of V 2 n calculated with a hard EOS have been multiplied by a factor of 0.58 (filled squares) 
to show the similarity with the shape of the SM predictions (open circles). One observes that the ratio 
of V 2 n versus time reflects closely the ratio of the integrated momentum changes due to the mean field 
gradients shown in Fig.4. 


Table 2: Passing time (see text) of the projectile and the target as a function of the incident energy in 
the +1®”^ Au collisions 


Einc. (A.MeV) 150 250 400 600 800 1000 1200 1500 

tpass (fm/c) 48.3 37.4 29.6 24.1 20.9 18.7 17.1 15.3 


scaling factor of 0.58 to the data for the hard EOS in Fig. [5l Despite the fact that the 
absolnte size of the signal differs by almost a factor of 2 for the EOS nnder considerations, 
this difference becomes stable fairly qnickly in the conrse of the reaction. The elliptic flow 
signals Vn themselves need rather long times to obtain their final (observed) shape. This 
is illnstrated in more detail in Fig. [6l the valne of V 2 at mid-rapidity stabilizes qnite fast, 
already at the fnll passing time, bnt it needs nearly twice as long to reach its final shape 
close to the target and projectile rapidities. 

We demonstrated that the mid-rapidity signal of V 2 {yo) develops rather rapidly within 
the IQMD model, stabilizing already after the passing time. This is an indicator that 
the spectator clock is working as anticipated. The elliptic flow at target and projectile 
rapidities takes a longer time to evolve since it originates also from late interactions with 
the spectator material. 

A decomposition of the elliptic flow into potential and collisional contribntions will 
elncidate the nnderlying reaction mechanisms. In Fig. [3 it is illnstrated how the hnal 
V 2 {yo) is obtained by NN collisions only, the mean held only, and by both operating 
together. The IQMD model allows to record the moment change due to collisions and due 
to potential separately for each time step. Integrating all contributions over the reaction 
time will yield the total momentum transfer vector due to the potential and analogously 
caused by the collisions. Both quantities have been analyzed to obtain ^ 2 ( 2 / 0 )- Note, that 


10 






0.1 


0.05 


-1 


1.5A.GeV, SM-t/t 

-0.39 

.0.78 

1.05 

. 1.63 

-3.92 


-0.5 0 


0.5 


y (u >0.4) 

■'O MO ' 


Figure 6: Results of IQMD calculations for the elliptic flow V 2 of protons as a function of scaled ra¬ 
pidity j/o at various values of the scaled time tjtpass (different line styles) and for an incident energy 
Ebeam =1.5 A.GeV. The protons are selected to have high transverse momenta Uto > 0.4 (see Section[2]). 
The IQMD calculations are the same as in Fig. |4] for a SM EOS. 


the values of those different ^ 2 ( 1 / 0 ) contributions cannot be simply added, since they do 
not reflect an absolute magnitude and the hnal momentum is obtained by a vector sum of 
the different contributions. 

We hnd the following features: 

a) The collisional contribution to V 2 alone is small and little affected by different EOS. 
This suggests that moderate medium changes of the NN cross sections that go beyond 
’zero order’ Pauli blocking corrections would not critically affect our conclusions on the 
EOS. Furthermore, the collisions by themselves cannot predict the measured elliptic flow. 
This illustrates the inadequacy of the cascade approaches for flow observables at these 
beam energies. 

b) The elliptic flow signal V 2 is more affected by the potential effects. One observes 
signihcantly differences of the ^ 2 ( 1 / 0 ) when applying different equations of states. 

Following this observation - that the elliptic flow signal is dominated by the acting 
potentials - we propose to construct a characteristic density which is connected to the 
elliptic flow of protons by calculating the mean value of the density seen by all protons 
with a final Uto > 0.4 weighted by the strength of the force due to the local mean field. 
This quantity integrated up to the full passing time is displayed by the line in the right 
hand-side of Fig. |8]as a function of bombarding energy. This variable is shown for the SM 
EOS which has been qualihed by the FOPI measurements. The shaded area represents the 
variances of the ” force-weighted” densities (corresponding to finite size effects). Using such 
a method we are not limiting ourselves to some small most central volume of the system a 
procedure which then naturally yields signihcantly larger densities. Those are depicted in 
the left hand-side of Fig. [HI The maximum density in the central volume max < p/po > {t) 
is shown as function of bombarding energy, as well as the maximum density seen by all 
protons with hnal ||/o| < 0.8 and the maximum density of protons with |?/o| < 0.8 and 
Uto > 0.4. In the innermost center of the collision densities up to 3.5po are reached at the 
highest beam energies measured by the FOPI collaboration. The density experienced by 
all protons measured in the region around mid-rapidity ||/o| < 0.8 is substantially lower 
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Figure 7: Elliptic flow V 2 of protons having finally a high transverse momentum uto > 0.4 (see Section [^l 
as a function of their scaled rapidity yo predicted by IQMD. Results for SM and HM EOS’s are shown 
at 0.6, 1 and 1.5A GeV incident energies at the late time t = 60 fm/c independent of energy. The total 
elliptic flow (’air, black lines) is presented along with its decomposition into the two contributions, i.e. 
mean held (’force’, dashed lines) and collisions (dotted lines). 


(up to 2po),and this value is not changed by the additional experimental cut Uto > 0.4. 

Our method of averaging with the strength of the force is in the same spirit as proposed 
in for the GMR. As can be seen in Fig. [HI the gap to the 2/3po data of the GMR 
could be closed by lowering the incident energies at least down to 0.1 A.GeV. Even if 
the densities extracted here are smaller than the maximum density reached during the 
collisions, the presence of the highly compressed matter (up to 3po) is mandatory for 
building up the high density gradients leading to the repulsion of particles and finally the 
elliptic flow signals. 


4. Comparison to microscopic calculations 


In view of the large number of proposed nuclear matter EOS in the literature, see for 
example the recent discussion of EOS derived using 240 (!) ’phenomenological’ (Skyrme 


interactions) approaches in 3^, we cannot here compare our constrained ’experimental’ 


EOS in Fig. 3 with all these various predicted behaviors, but hope that future theoretical 
efforts will consider our results for the density range up to p/po = 2.5. 

By ’microscopic’ calculations we define here theoretical efforts to calculate the inhnite 
nuclear matter EOS using nucleon-nucleon and pion-nucleon scattering data and some 
information on multi-body forces (3N, 4N) from light nuclear clusters, including of course 
the deuteron, without introducing further ad hoc parameters. Taranto et ah have 
recently tested five such models for their consistency with known astrophysics and nuclear 
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Figure 8: Left side: The maximum of the density < p(t)/p >max reached in the central volume of the 
collision (shown in Fig. 0] for 1.5A GeV) as a function of beam energy is shown as open triangles, the 
maximum of the density < p{t)/pQ >y-cut for protons finally in the rapidity range |j/o| < 0.8| as reversed 
triangles and as open circles the maximum of < p{t)/po >y-ut-cut seen by protons with |j/o| < 0.8| and 
Uto > 0.4. Right side: ”Force-weighted” mean value of the reduced density (see text) seen by protons 
Uto > 0.4 and |j/o| < 0.8 in their final configuration averaged until the passing time. The error bars 
represent the standard deviations of the distributions. 


physics constraints. They came up with open questions concerning the predictions for high 
densities (p > 3po) that one needs to know in compact star modeling. An intriguing point 
was that three of the models violated causality at the predicted central densities and the 
suggestion that even Dirac-Brueckner-Hartree-Fock (DBHF) calculations may be in error 
due to approximations made in order to make the calculation tractable. 

We choose here to confront three representative microscopic calculations with our new 
constraints. First a DBHF calculation of both symmetric and pure neutron matter 
which was presented in a conveniently parameterized form allowing to estimate with use 
of the parabolic assumption (involving S'^) the deviation of the EOS in going from strictly 
symmetric to Au-like = 0.0391) composition, where 6 = {N — Z)/{N + Z). The second 


comparison uses work from the ’Idaho group’ [37| allowing to compare the historically 
earlier DBHF approach with the more recent chi rally motivated theories. Finally, we 
compare with two different EOS published in 2005 |38| one having the interesting feature 
of including the possibility of (virtual) A excitations explicitly. 

The EOS calculated with DBHF in using the Bonn A jsll nucleon-nucleon potential 
(full green curve) is confronted in Fig. [9] with our constrained FOPI-IQMD EOS (shaded 


band around the median black curve). The authors 3^ have made efforts to check the 


validity of averaging approximations often used in DBHF to reduce the computational 
cost. 

As mentioned before we have also estimated the ’correction’ to the symmetric EOS due 
to the degree of isospin asymmetry in gold (An), see the dashed curve. It is located roughly 
1 — 1.5A MeV above the symmetric system. It is visible that the shape of the ’An’ EOS is 
not dramatically different from the strictly symmetric matter EOS. Of course the heavy ion 
flow data, being sensitive only to the pressure, i.e. essentially the derivative of such curves. 
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Figure 9: Symmetric nuclear matter EOS from ref. 36| using Bonn A potential (full green curve). The 
dashed red curve is the EOS for Au-like systems using neutron matter calculations from [s^ . The shaded 
error band surrounding the black curve is the FOPI-IQMD constraint. 


does not allow to determine the absolute vertical position of the energy/nucleon EOS. The 
uncertainty of this position is roughly indicated by the rectangle in the figure. As there 
seems no clear consensus on the widths of this ’saturation rectangle’ in the literature we 
have adopted here approximately the sizes given in two theoretical papers, If one 


adopts conclusions from given observables (such as nuclear masses and electron-scattering 
data) one can come to tighter, but still model dependent limits. As just one example, 
one obtains very tight limits (<0.2 MeV) on the bulk binding energy if one interprets the 


mass-number proportional coefficient as holding for inhnite NM, see a recent effort in ^2 


The narrowing ’error’ band of the empirical curve in approaching the saturation point 
is a consequence of constraining the curve to cross the most probable saturation value and 
ignores the shown uncertainty of this very fundamental point. 

Constraining the stiffness of pure neutron matter at densities significantly exceeding 
the saturation density is and will be a difficult task as realistic heavy ion systems, even 
at radioactive beam facilities, will not exceed 6^ values larger than 5%. It is therefore not 
surprising that information relevant to compact star physics ji^ is so far primarily limited 
to subsaturation densities that are accessible by other means. 

In Fig. [10] we show two symmetric nuclear matter EOS from ref. j^. One, full 
triangles, labeled ’DBHF’ in the hgure, is using a meson theoretic potential together with 
the DBHF method, the other, open triangles, labeled ’Chiral’, makes use of effective held 
theory (EFT) with density dependent interactions derived from leading order chiral three- 
nucleon forces. 

One can see that the DBHF curve tends to become stiffer at densities exceeding p > 
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Figure 10: Symmetric nuclear matter EOS after ref. [37|. Full (red) triangles: DBHF calculation, open 


(blue) triangles: chiral effective field theory. The shaded band is the FOPI-IQMD constraint. 


1.5po- On the other hand the ’chiral’ version seems to match close to perfectly with the 
FOPI-IQMD constraint (shaded). As the EFT approach, respecting QCD symmetries, is 
a low-momentnm theory sorting the contributions of various order according to powers of 
momentum the question arises up to what density one can trust the predictions. In 
Machleidt pointed out that chiral NN potentials do not make any reasonable predictions 
beyond 300 MeV laboratory energy. He ends up concluding that ’one may trust the chiral 
perturbation theory up to densities around 4po’- 

The treatment of the short range parts of the nucleon-nucleon forces outside of QCD 
lattice theory is still subject to different approaches, not only between the above mentioned 
DBHF and Chiral EOS, but also in some other ’chiral approaches to nuclear matter’ (j^ . 
As we wish to limit the scope and size of the present contribution, we will not describe or 
comment these various approaches. However we can show, see Fig. [TT|, how the theoretical 
investigations may use the FOPI-IQMD constraint as rough guide line. 

Using the chiral approach, the authors predict two rather different EOS depending 
on whether they include or not virtual A excitations. See the hgure caption. The predic¬ 
tions including the A are in agreement with the empirical saturation point. The theory is 
not completely free of some adjustments that go beyond the use of nN and NN scattering 
data. The short-range parameters of the theory have been ’fine-tuned’ to reproduce the 
(—16) MeV binding energy per nucleon at saturation. However, the virtual A-excitations 
then help to locate the EOS at the right horizontal place around p = 0.16 fm“^. A look 
at Fig. [TT] shows that the A leads to a rather marked stiffening of the EOS {Kq = 304 
MeV) projecting the curve outside the shaded FOPI-IQMD constraint. It is not obvious 
to us what this discrepancy means. It could be that the apparent ’heavy-ion EOS’ actually 
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Figure 11: Symmetric nuclear matter EOS after ref. [s^. The dashed (blue) EOS curve includes pions 
and nucleons as active degrees of freedom. The full (red) curve includes effects from 2t: exchange with 
virtual A-excitation. The shaded band is the FOPl-lQMD constraint. 

differs from a ’cold’ EOS because, due to the finite temperature in the reaction the A are 
real rather than virtual. The theoretical ’A stiffness’ could then be a dispersion effect 
rapidly changing with temperature. If so, it would be eminently important to pursue this 
question. Naively, one would expect that the (real) possibility to store energy internally 
into the nucleon (at high T) would lower the kinetic energies and hence the pressure and 
stiffness, some kind of mild ’phase transition’ when new degrees of freedom open. 

5. Discussion 

We have shown in section [2]that a single parameter which we dubbed n 2 n, characterizing 
elliptic flow over a large rapidity interval (80% of the initial gap) allows to discriminate 
rather clearly between a soft (SM) and a stiff (HM) EOS. V 2 n differs by a factor 1.63 
over an energy range extending from OAA GeV to 1.5A GeV when using IQMD, while 
the experiment allows to constrain V 2 n within roughly a factor 1.1. The relevant density 
range was estimated from the simulations to span p^ax = (1-1 — 3.0)^0, where the high 
densities reached during the collisions are important for the evolution of the hnal elliptic 
flow pattern. This makes the ’flow method’ competitive and complementary with the ’kaon 
method’: the latter reaches such a high SM/HM discrimination only around 0.8A GeV far 
below threshold. It is gratifying that both methods lead to the same conclusion, however. 

In Fig. [12] the range of pressure and density assessed in this work is shown by the 
area enclosed by the dashed lines. The pressure for symmetric matter at zero temperature 
predicted by the FOPI-IQMD constrained EOS with k, = 190 ± 30MeV is shown as 
magenta line. To put these results in context with others they are compared to the 
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Figure 12: Nuclear EOS in terms of pressure versus density. The pressure as deduced from the soft EOS 
which was used in the IQMD model calculations is shown by the magenta line. The dashed curves enclose 
the error band in the determination of the relevant pressure region. The area framed with lines originates 
from the analysis of Q 


constraints reported in the analysis of [8[ at higher energies and respective densities. Onr 
data is in agreement with the earlier hndings giving slightly more restrictive constraints 
towards a more soft nnclear matter EOS at higher densities. 

Convincing conclnsions on basic nnclear properties reqnire a snccessfnl simnlation of 
the fnll set of experimental observables with the same code nsing the same physical and 
technical parameters. For a nnmber of observables this goal has been reached, for some 
other data this is not yet the case. 

It is well known that in the beam axis system elliptic flow is not ind^endent of the flow 
axis and hence the magnitude of the directed flow. In contrast to ref. |8| , the IQMD model 
reproduces the directed flow using the same EOS parameters. In ref. this is shown 
(Fig. 49) for directed flow of deuterons in the same beam energy and centrality range. 
Here too the experimental errors were signihcantly smaller than the HM/SM difference. 

Furthermore, in the same work. Fig. 51, radial flow of the light clusters was well 
reproduced (and found to be insensitive to EOS). 

Stopping information 47], l48| and Fig. 52 in 22| relevant to the question of in-medium 
nucleon-nucleon cross sections, requires still some hne tuning on the 20% level. Incomplete 
stopping could mockup a soft EOS if it is incorrectly accounted for in the simulations (45l |. 

Pion yields, measured with 10% accuracy, ji^, differ only by about 10% between the 
HM and SM options, and the yields were found to have about 20% lower values than the 
SM option of IQMD. It should be kept in mind that the elementary pion cross sections 
are known with limited precision. 


17 










Some less desirable features of our version of IQMD will not be hidden: 

a) the momentum dependencies of protonic vi are not correctly reproduced, while 
it is increasingly correct the heavier the cluster; this could be due to a lack of quantum 
fluctuations in a semiclassical code; we expect that fluctuations have a diminished influence 
on V 2 n because it is determined by a large part of phase space; 

b) the degree of clusterization is underestimated; this deserves further investigations. 

It is worthwhile mentioning that in some comparisons were also done with lower 

energy data. At 0.15A GeV the difference HM/SM was found to be rather small. The 
deuteron elliptic flow was well reproduced by IQMD, Fig. 38. The small sensitivity to 
EOS at beam energies well below OAA GeV is not surprising in view of the flatness of the 
energy/nucleon EOS (see Fig. 3) in this regime, implying relatively low pressures. We 
also remark that at 0.25A GeV the Vi flow of all light clusters was well reproduced. The 
reproduction of low energy data, where the effect of the EOS is still moderate, is useful to 
ensure the adequacy of some technical parameters of the code. 

Last, but not least, essentially the same IQMD code (except for some kaon specific 
sections) was used for the ’kaon method’. 

Despite these encouraging results, there are two important points to be aware of: 

1. data that are claimed to constrain such fundamental properties as is the EOS of 
nuclear matter require conhrmation by independent experimental efforts 

2 . there is an urgent need to prove that the conclusions using a particular transport 
code, here IQMD, are not limited to using this particular code. 

We shall briefly comment on these two points. 

Before the ’experimental’ EOS is finally ready for textbooks one needs a general con¬ 
sensus on the measured data with use of various apparatus. At the time of writing this 
goal is not yet reached in a satisfactory way. In ref. Q] two of the at that time available 
experimental data points on elliptic flow within the SIS/BEVALAG energy range (below 
2A GeV) survived the critical eye of the authors: a point (|n 2 | ~ 7.7%) around 0.4A GeV 
from the Plastic Ball Gollaboration ^ and a point near I.IA GeV (6.2%) from the EOS 
Gollaboration both obtained at the BEVALAG accelerator in Berkeley (USA) with a dif¬ 
ferent apparatus. We could not hnd the (non-rotated) Plastic Ball value in the literature 
as the authors of refs. choose at the time the so called ’flow axis’ as reference to 

determine elliptic flow. To our knowledge the AGS data were not ’rotated’ (see 0 
for a brief discussion of the issue of rotated reference systems). A close comparison in ^ 
with the Plastic Ball’s rotated data suggested disagreement with the FOPI data in at least 
part of the overlapping incident energy range. The 1.1 A GeV point accepted in was 
hrst cited in , but we could not hnd a detailed documentation in the refereed literature 
on how possible apparatus cuts and other effects were taken into account. Optimistically, 
the l.lA GeV value is compatible with the FOPI data if one assumes a compensation of 



The second important point to be discussed is that the conclusions drawn from the 
comparison with simulation data obtained from various transport codes must be consis¬ 
tent. There are several issues that need further efforts by the community: momentum 
dependencies (including clean Lorentz covariance at beam energies exceeding lA GeV), 
clusterization and entropy balances, in-medium nucleon-nucleon reactions etc. Efforts to 
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bring transport code experts together, to see if a convergence of analyses results can be 
reached, have been done but urgently need to be continued including energies up to 
FAIR accelerator design. 

What can we say about the future besides asking for efforts to reach consensus on both 
data and transport code analyses? 

The sensitivity of V 2 to the assumed stiffness of the EOS was shown to diminish 
rapidly below OAA Gev and above 6A GeV. We now tend to associate the weak sensitivity 
at low energies to the combined action of insufficient stopping ^ and still relatively 


modest compression. At energies well above 4AGeV, on the other hand, the shadowing 
effect by spectator nucleons (and therefore the time clock) is no longer very effective as 
suggested by the observed sign change of V 2 at AA GeV [Ml]- Presumably the spectator 
clock can be used to extend improved EOS constraints to densities (3-4 po) higher than 
was possible in the present work (limited to beams of 1.5A GeV incident energy) in future 
accelerator systems such as FAIR in Darmstadt, Germany. However it is clear that beyond 
4A GeV, other ideas are needed to extract EOS information from heavy ion data. 

As a main conclusion of the present work we believe we can say that the feasibil¬ 
ity of establishing reasonably tight empirical constraints on the nuclear EOS has been 
demonstrated. 
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